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<^ ! ABSTRACT 

■ 

' We present the first results of the unbiased survey of the L1157-B 1 bow shock, obtained with HIFI in the framework of the key program Chemical 
(-H Herschel surveys of star forming regions (CHESS). The L1157 outflow is driven by a low-mass Class protostar and is considered the prototype 
^ ^ of the so-called chemically active outflows. The bright blue-shifted bow shock B 1 is the ideal laboratory for studying the link between the hot (~ 
I 1000-2000 K) component traced by H2 IR-emission and the cold (~ 10-20 K) swept-up material. The main aim is to trace the warm gas chemically 
O enriched by the passage of a shock and to infer the excitation conditions in LI 157-Bl. A total of 27 lines are identified in the 555-636 GHz region, 
^ down to an average 3cr level of 30 mK. The emission is dominated by CO(5-4) and H20(lio-loi) transitions, as discussed by Lefloch et al. in 
C/3 , this volume. Here we report on the identification of lines from NH3, H2CO, CH3OH, CS, HCN, and HCO^. The comparison between the profiles 
Ci , produced by molecules released from dust mantles (NH3, H2CO, CH3OH) and that of H2O is consistent with a scenario in which water is also 
formed in the gas-phase in high-temperature regions where sputtering or grain-grain collisions are not efficient. The high excitation range of the 
^ observed tracers allows us to infer, for the first time for these species, the existence of a warm (> 200 K) gas component coexisting in the B 1 bow 
^ structure with the cold and hot gas detected from ground. 
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• 1. Introduction molecular bow-shock is essential bevause it represents a very 

powerful diagnostic tool for probing its physical conditions. 
' A newborn protostar generates a fast and well coUimated let, t.. t 1 1 .a 1 . j . j- . . j . ^. ^. 

.,, J ji -J 1 • J T . • 1 J The LI 157 outflow, located at a distance estimated to be be- 

possibly surrounded by a wider angle wind. In turn, the ejected /, , i^rvrv-rt j /ir • nrrrrn. 

^ ■ , . • /, X . 1 . 11- ... J- tween 250 pc (Looney et al. I2007I I and 440 pc (Viotti 119691 ) 

^-H material drives (bow-)shocks travelling through the surrounding , . . . • j . ,i — r . • , r 

. J- J J 1 fi -1 ■ 1 .■ may be regarded as the ideal laboratory tor observing the er- 

high-density medium and traced by H2 ro-vibrational lines at „ „ , , , , . , • , , r , 

. / ^ J. J r,nnA 1^ A rects or shocks on the gas chemistry, being the archetype or the 
excitation temperatures or around 2000 K. As a consequence, .. j . ■ „ a /V. . -^i . ■<rck^ 



X 



1 J u /irv ^r. 1 1 a r ji so-callcd chcmically Hch outflows (Bachillcr ct al. 200 1 ). Thc 

slower and cold (10-20 K) molecular outflows are formed by , , . ,a i t \ r^ . . m a oarvoo^ ^^ci 

,^ . , . , \ „„ c^u , • ^ ^i. low-luminosity (4-11 Lo) Class protostar IRAS20386+6751 

H swept-up material, usually traced by CO. Shocks heat the gas , . . ^ . . . • j • . 

rrt j^- 1 i_ j.u u i drives a precessing powerful molecular outflow associated with 

I and trigger several processes such as endothermic chemical re- , , ^ . , . nr^rm j- m .t 

,■ . 1 1 .■ ■ • n 1 several bow shocks seen in CO (Gueth et al.'199o^ and in IR Ht 

actions and ice grain mantle sublimation or sputtering. Several . /r^ • o r-- , -rc .rniTTji xt r u . ^r.r,r. t 

, , . , ■ -r- , . , . images (Davis & Eisloffel l 19951 Neufeld et al. 2009). In particu- 

molecular species undergo signmcant enhancements in their , , , . , , . ' , , , , „ .m /r-- m . 

, , K u 1 o ni 1 i .nnot i. j lar, the brightest blue-shifted bow-shock. Called B 1 (Fig. lib, has 

abundances (see e.g.. van Dishoeck & BIake |1998l ). as observed ' . , j ■, , run j ^^T a • r- 

, , . , , 7 , r been extensively mapped with the PdB and VLA interferome- 

by observations at millimeter wavelengths towards a number of j . • .■ • . j . 

r, , ,rvr,o Ti 1 •.. D rr ^ , r.r.-, tei"s at mm- and cm-observations revealing a rich and clumpy 

outflows (Garay et al. 1 998; BachiUer & Perez Gutierrez 1997, , , , • , . . „ i- . • 

TiTinTu ' 1 1 r^rvrv-i. t-u v i i. . structurc, the clumps being located at the wall of the cavity 

BP97 hereafter; j0rgensenetaLKOO7). The link between the gas ., , , ^ „ o „ un i mng i .u . i inr,^ 

i,?T^ J . . i j With an arch-shape (Tafalla & Bachiner |1995l Gueth et al. 19981 

components at ~ 10 K and the hot 2000 K shocked component „ j „• • . i i-^Am i u n-^rr^nn j n . i ^nrv<^ 

^ . , , J- 1 , ^^ ■ J r Bcnedettim et al. 120071 hereafter BVC07; Codella et al. 20091). 

is crucial to understanding how the protostellar wind transfers , , ic-, m • ..T j . . . . . . . j . 

, 1. 1 . T.- . T .u- LI 157-Bl IS well traced by molecules thought to be released by 
momentum and energy back to the ambient medium. In this con- , , , ,, ^tt ^rr j xm .. • 
, , ..-^ ^ 1 . • 1 • • r • 1 dust mantles such as HiCO, CH3OH, and NH3 as well as typi- 
text, the understanding of the chemical composition of a typical , r, ■ , ji i i n-^/ ^ji^ . 
° ^ -^^ cal tracers of high-speed shocks such as SiO (e.g., Gusdorf et al. 

* Herschel is an ESA space observatory with science instruments 2008 ). Temperatures ^ 60-200 K (from NH3, CH3CN, and SiO) 

provided by European-led principal Investigator consortia and with im- as well as around 1000 K (from H2) have been de rived ( Tafalla & 

portant partecipation from NASA Bachiller il995i Codella et al. i2009i Nisini et al. 120071 in prep.). 
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Fig.l. The Bl clump. PdBI emission of CH30H(2i-li)A" 
(grey) on the CS(2-1) one (contours), from BVC07. The maps 
are centred on the coordinates used for the present HIFI observa- 
tions ajzooo = 20'^ 39™ l(f.2, (Jj2ooo = +68° 01' 1T5, i.e. at Aa = 
+25'.'6 and A5 = -63'.'5 from the driving protostar. The labels in- 
dicate the main Bl clumps detected in different tracers. Circles 
are for the HPBWs of the HIFI data presented here (39") and 
of Band 7 (11"), i.e., at the highest frequencies of the CHESS 
surveys. 

However, a detailed study of the excitation conditions of the B 1 
structure has yet to be completed because of the limited range of 
excitation covered by the observations performed so far at cm- 
and mm- wavelengths. Observations of sub-mm lines with high 
excitation (> 50-100 K above the ground state) are thus required. 

As part of the Herschel Key Program CHES^l (Chemical 
Herschel Surveys of Star forming regions), L1157-B1 is cur- 
rently being investigated with an unbiased spectral survey us- 
ing the HIFI instrument (de Graauw et al. 2010). In this Letter, 
we report the first results based on HIFI observations in the 555- 
636 GHz spectral window, confirming the chemical richness and 
revealing different molecular components at different excitation 
conditions coexisting in the Bl bow structure. 

2. Observations 

The observations were performed on 2009, August 1, during the 
Performance Verification phase of the HIFI heterodyne instru- 
ment (de Graauw et al. 2010) on board of the Herschel Space 
Observatory (Pilbratt et al. 2010). The band called lb (555.4- 
636.2 GHz) was covered in double-sideband (DSB) with a total 
integration time of 140 minutes. The Wide Band Spectrometer 
was used with a frequency resolution of 1 MHz. The typical 
HPBW is 39". The data were processed with the ESA-supported 
package HIPlfl (Herschel Interactive Processing Environment) 
for baseline subtraction and sideband deconvolution and then 
analysed with the GILDA^ software. All the spectra (here in 
units of antenna T^) were smoothed to a velocity resolution of 
1 km s"\ except those showing the weakest emission, which 
were smoothed to lower spectral resolutions (up to 4 km s"'). 
At a velocity resolution of 1 km s the rms noise is 6-13 mK 

' http://www-laog.obs.ujf-grenoble.fr/heberges/chess/ 

^ HIPE is a joint development by the Herschel Science Ground 

Segment Consortium, consisting of ESA, the NASA Herschel Science 

Center, and the HIFI, PACS and SPIRE consortia. 
^ http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 3. Top and Middle panel: Comparison between the profiles 
of NH3(lo-Oo), multiplied by a factor 7.4, H2CO(8i7-7i6), mul- 
tipled by a factor 22.5, HCN(7-6), multipled by a f actor 9.0, and 
H20(lio-loi), the latter from Lefloch et al. (120 1 01) . The vertical 
solid line indicates the ambient LSR velocity (+2.6 km s^'). The 
velocity ranges arbitrarily defined as HV (-20,-6 km s"'; traced 
by H2O), MV (-6, -1 .5 km s ' ; outlined by the HCN and H2CO 
secondary peak), and LV (-1.5,-1-2.6 km s"'; the rest of the blue 
wing) are drawn (see text). Bottom panel: Intensity NH3/H2O 
line ratio as a function of velocity. 



(Ta scale), depending on the line frequency. The main-beam effi- 
ciency (ri„b) has not yet been reliably determined. When needed, 
we adopted an average rj^b of 0.72. 

3. Different tracers at different velocities 

A total of 27 emission lines were detected, with a wide range 
of upper level energies, from a few tens to a few hundreds of 
Kelvin. Table 1 lists the spectroscopic and observational param- 
eters of all the transitions. For the first time, high excitation 
(up to ^ 200 K) emission lines related to species whose abun- 
dance is largely enhanced in shocked regions were detected. The 
CO(5^) and H20(lio-loi) lines are analysed in Lefloch et al. 
(2010). Figure 12] presents representative examples of line pro- 
files observed towards L1157-B1. All the spectra contain lines 
with blue-shifted wings peaking near km s \ which have a 
terminal velocity equal to ~ -8,-6 km s Previous PdBI ob- 
servations showed that L1157-B1 is associated with very high 
velocities (HVs) of as low as ^ -20 km s"' (v'lsr - +2.6 km 
s BP97). We cannot exclude the lack of detected emission in 
the HV regime in the present HIFI spectra being caused by their 
relatively low signal-to-noise (S/N) ratio. The PdBI images in- 
dicate that the brightness of the emission lines in the HV regime 
is indeed weaker than the emission at low velocities by a fac- 
tor of 5-10. The spectra in Fig. |2] clearly show that this weak 
emission would lie below the noise. On the other hand, the HV 
gas is detected in the very bright lines of CO and H2O (Lefloch 
et al. 120 1 Oi l. We note that the HV emission is mostly confined 
to within the eastern Bla clump (Fig.[T]i, within an emitting re- 
gion of size < 10" (Gueth et al. IT^Sl BVC07), whereas low 
velocity lines originate in both the bow-structure and the walls 
of the outflow cavity (e.g., the BOe and BOd in Fig.[T]i, of typical 
size 15"-18". Therefore, the forthcoming HIFI-CHESS obser- 
vations at higher frequencies and higher spatial resolution (see 
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Table 1. List of molecular species and transitions observed with HIFI (Band lb): CO and H2O emission is discussed in Lefloch et 
al. (1201 Oi l. Peak velocity and intensity (in scale), integrated intensity (Fint), as well as the terminal velocities of the line emission 
(Vmin and Vmax) are reported. 



Transition 


vo 




^peak 


rms 


Vpcak 










(MHz) 


(K) 


(mK) 


(mK) 


(km s-') 


(km s-') 


(km s-') 


(Kkms-1) 


o-HjOd.o-loi) 


556936.002 


27 


910(17) 


17 


-0.37(1.00) 


-25.4 


+7.6 


11.68(0.10) 


CH,OHE(ll29 - IO19) 


558344.500 


168 


47(10) 


10 


+0.60(1.00) 


-2.9 


+2.7 


0.16(0.02) 


o-H2CO(8i8-7i7) 


561899.318 


118 


92(4)* 


8 


+0.38(0.14)* 


-4.5 


+3.6 


0.49(0.03)* 


CH30HE(32.2-2i,2) 


568566.054 


32 


42(8) 


8 


+0.60(1.00) 


-2.6 


+2.7 


0.31(0.02) 


o-NHjdo-Oo) 


572498.068 


28 


122(7) 


7 


+ 1.03(0.80) 


-6.9 


+5.0 


0.89(0.03) 


CO(5-4) 


576267.931 


83 


883(10) 


10 


+ 1.60(1.00) 


-36.5 


+6.0 


49.30(0.07) 


p-H2CO(8o8-7o7) 


576708.315 


125 


41(7) 


7 


+ 1.80(1.00) 


-5.4 


+2.7 


0.22(0.02) 


CH3OHA- (22,1 -l,,o) 


579084.700 


45 


53(8) 


8 


+0.60(1.00) 


-6.0 


+3.7 


0.29(0.03) 


CH30HE(12u2 - lli.ii) 


579151.003 


178 


44(8) 


8 


-0.30(1.00) 


-4.0 


+4.0 


0.25(0.02) 


CH30HA+(12o,,2-llo.ii) 


579459.639 


181 


48(7) 


7 


+0.90(1.00) 


-3.9 


+4.9 


0.27(0.02) 


CH3OHA+ (22.0-1 1.1) 


579921.342 


45 


42(7) 


7 


-0.60(1.00) 


-3.3 


+2.6 


0.21(0.02) 


CH30HE(122 10 - II29) 


580902.721 


195 


11(4) 


4 


-3.00(3.00) 


-3.9 


+2.5 


0.09(0.02) 


CH3OHA+ (61.6 -5o.5) 


584449.896 


63 


88(9) 


9 


-0.30(1.00) 


-6.0 


+2.9 


0.58(0.03) 


CS(12-11) 


587616.240 


183 


23(2)* 


5 


-0.61(0.57)* 


-7.3 


+6.5 


0.19(0.03)* 


CH3OHA+ (735-624) 


590277.688 


115 


42(9) 


9 


+0.60(1.00) 


-1.0 


+3.0 


0.16(0.02) 


CH3OHA- (734-625) 


590440.291 


115 


40(9) 


9 


-0.60(1.00) 


-2.6 


+0.8 


0.10(0.02) 


CH30HE(9o.9-8i,8) 


590790.957 


110 


40(10) 


10 


+0.60(1.00) 


-A.l 


+4.3 


0.28(0.03) 


o-H2CO(8i7-7i6) 


600374.604 


126 


29(7)* 


9 


-0.14(0.57)* 


-3.0 


+ 1.8 


0.19(0.04)* 


CH30HE(42 3 - 3i3) 


616979.984 


41 


47(6) 


6 


+0.90(1.00) 


-4.5 


+2.6 


0.26(0.02) 


HCN(7-6) 


620304.095 


119 


94(7) 


7 


-0.60(1.00) 


-7.6 


+3.2 


0.68(0.03) 


HCO+(7-6) 


624208.180 


119 


30(3)* 


8 


+0.53(0.47)* 


-3.6 


+4.5 


0.11(0.03)* 


CH3OHA- (32 2 - 2i 1) 


626626.302 


52 


18(5) 


5 


-1.20(4.00) 


-2.5 


+6.3 


0.18(0.07) 


CH30HE(13i.i3 - 121.12) 


627170.503 


209 


33(7) 


7 


-0.30(1.00) 


-3.2 


+2.5 


0.15(0.02) 


CH30HA+(13o.i3-12o.,2) 


627558.440 


211 


41(9) 


9 


+0.60(1.00) 


-3.6 


+2.5 


0.19(0.02) 


CH3OHA+ (3i2 -2,2) 


629140.493 


52 


52(9) 


9 


+0.60(1.00) 


-3.5 


+3.7 


0.24(0.03) 


CH3OHA+ (7,,7-6o,6) 


629921.337 


79 


70(13) 


13 


+ 1.50(1.00) 


-3.9 


+3.7 


0.37(0.04) 


o-H,CO(9i9-8i8) 


631702.813 


149 


66(4)* 


9 


+0.45(0.17)* 


-1.5 


+2.6 


0.22(0.02)* 



" Frequencies and spectroscopic parameters have been extracted from the Jet Propulsion Laboratory molecular database (Pickett et al.[T998j for 
all the transition except those of CH3OH, which have been extracted from the Cologne Database for Molecular Spectroscopy (Miiller et al. 12005b . 

Upper level energies refer to the ground state of each symmetry. * Gaussian fit. 
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Fig. 2. Molecular line profiles observed towards LI 157-Bl: species and transitions are reported in the panels. The vertical solid line 
indicates the ambient LSR velocity (+2.6 km s ' from C''^0 emission; BP97), while the dashed one is for the secondary peak at 
-4.0kms-i. 



the dashed circle in Fig. [TJ should allow us to study the HV 
wings in species other than CO and H2O. 

The uniqueness of HIFI lies in its high spectral profile reso- 
lution for many high excitation transitions of a large number of 
molecular species. The analysis of the present HIFI spectra re- 
veals a secondary peak occuring between -3.0 and ^.0 km s"' 



(here defined medium velocity, MV) and well outlined by e.g., 
HCN(7-6). The MV peak is also visible in NH3(lo-0()) and in 
some lines of CH3OH and H2CO (see Fig.|3]l, but its occurrence 
does not show any clear trend with the choice of tracer of line 
excitation. No single-dish spectra had previously detected this 
spectral feature (BP97; Bachiller et al. l2001l l. An inspection of 
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the spectra observed at PdBI shows that the MV secondary peak 
is observed in a couple of lines of the CH30H(2k-1k) series (see 
Fig. 3 of BVC07) and only towards the western Bib clump (size 
~ 5")- This finding implies that there is a velocity component 
originating mainly in the western side of Bl, while the HV gas 
is emitted from the eastern one (see above). 

Figure [3] compares the profiles of the NH3(lo-Oo) and 
H2CO(8i7-7i6) lines with the H20(lio-loi) profile, where the 
S/N allows such an analysis (MV and LV ranges). By assuming 
that the emission in the MV range is optically thin (including the 
H2O line) and originates in the same region, we obtained from 
the comparison of their profiles a straightforward estimate of the 
relative abundance ratios of the gas at different velocities. As 
a notable example, the NH3/H2O intensity ratio decreases by a 
factor ~ 5 moving towards higher velocities (Fig. O, implying 
that a similar decrease in the abundance ratios occurs. This may 
reflect different pre-shock ice compositions in the gas emitting 
the MV emission. Alternatively, this behavior is consistent with 
NH3 being released by grain mantles, but water both being re- 
leased by grain mantles and, in addition, copiously forming in 
the warm shocked gas from endothermic reactions, which con- 
vert all gaseous atomic oxygen into water (Kaufman & Neufeld 
1997; Jimenez-Serra et al. 2008. and references therein). The 
water abundance may be enhanced with respect to ammonia in 
the fast and warm (> 220 K) gas, which might explain why the 
H2O wings are larger than those of NH3, CH3OH, and H2CO, 
all species being directly evaporated from dust grain mantles. 

4. Physical properties along the B1 bow shock 

We detected several lines from CH3OH (17 lines with upper 
level energies up to 211 K). We can derive a first estimate of 
the emitting gas temperature by means of the standard analysis 
of the rotational diagram. We show the case of methanol (A- and 
E-forms) in Fig. |4] The derived rotational temperature (Trot) is 
106 K (with an uncertainty of ~ 20 K), which represents a lower 
limit to the kinetic temperature (r^n). In the same figure, we re- 
port the methanol lines (2k-1k) observed with PdBI and whose 
intensity is integrated in the HIFI 39" beam. The T^ot derived 
from the ground-based data (based only on lines with < 50 
K; BVC07) is definitely lower, ~ 12 K, in perfect agreement with 
that found with the 30-m spectra in the same excitation range by 
Bachiller et al. (1995). As discussed by Goldsmith & Langer 
(II999I 1. this behavior may be caused by either two components 
at different temperatures or both non-LTE effects and line opac- 
ity. These two possibilities cannot be distinguished based only 
on the rotational diagram. However, given that a range of Tkin 
and is naturally expected in a shock, if we were to assume 
that two gas components provide an explanation, they would not 
only have different temperatures but also a different column den- 
sities. Taking the filling factor ff = 0.13, derived by the CH3OH 
maps obtained at the PdBI, the low temperature component col- 
umn density is 8 x 10'"* cm"^ (in agreement with Bachiller et al. 
11995) . whereas the high temperature component has a column 
density of around 10'"* cm"^. We note that the rotation diagrams 
obtained for the MV and LV CH3OH emission separately do not 
allow us to infer any clear difference. 

It is possible to more tightly constrain the emitting gas tem- 
perature and «Ht density for the species where the coUisional rate 
coefficients are known, by performing of a non-LTE analysis. To 
this end, we used the non-LTE excitation code RADEX with an 
escape probability formalism for the radiative transfer (Van der 
Ta k et al . 120071) coupled with the LAMDA database (Schoier et 
al. I2OO5I 1. Methanol is the species detected in the largest num- 
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Fig. 4. Rotation diagrams for the CH3OH transitions measured 
with HIFI (triangles) and from ground (PdBI; squares). Black 
and blue points are for A- and E-form, respectively. The param- 
eters A^u, and are, respectively, the column density, the 
degeneracy and the energy (with respect to the ground state of 
each symmetry) of the upper level. The derived values of the ro- 
tational temperature are reported: (i) 106 K, for the HIFI lines 
covering the = 32-211 K excitation range and (ii) 12 K (as 
Bachiller et al. ll995l l. for the PdBI lines, at lower excitation. 



ber of lines. The full non-LTE study will be reported in a forth- 
coming paper. Here we analysed only the E-form, for which the 
collisional rate coefficients are available (Pottage et al. '20041I. 
The major result of this analysis is that for a range of densities 
of lO-'-lO^ cm"-', the gas temperature exceeds 200 K. A similar 
result is obtained by considering H2CO emission. 

Finafly, by combining the HIFI CS(12-1 1) fine with CS(2- 

I) and (3-2) lines observed with ground-based telescopes, we 
also derive a kinetic temperature that is definitely above 300 
K for the outflowing gas. In this case, caution should be taken 
since we are abl eto trace different gas components, as suggested 
by CH3OH, the gas at higher excitation being traced by CS(12- 

II) . If we analyse only the (2-l)/(3-2) intensity ratio, the non- 
LTE approach does not allow us to constrain the temperature in 
this way, but we are able to infer «h, of around 4 x 10"^ cm"^. 
Interestingly, when we check for a possible dependence of 

on velocity, the LV range is found to be indicative of a denser 
medium (~ 10^ cm""*) by an order of magnitude with respect to 
the MV gas. 



5. Conclusions 

We have presented the HIFI unbiased spectral survey in the 555- 
636 GHz band towards the bright bow-shock Bl of the LI 157 
protostellar outflow. For the first time, we have detected high- 
excitation (up to ^ 200 K) emission lines of species whose abun- 
dance is largely enhanced in shocked regions (e.g., H2O, NH3, 
H2CO, CH3OH). This has allowed us to trace with these species 
the existence of a high excitation component with Tkin > 200- 
300 K. Temperature components from ~ 300 K to ~ 1400 K have 
been inferred from the analysis of the H2 pure rotational lines 
(Nisini et al., in prep.). Therefore the present observations pro- 
vide a link between the gas at r^n 60-200 K previously observed 
from the ground and the warmer gas probed by the H2 lines. We 
plan to perform additional HIFI observations in the THz region 
towards LI 157-Bl to observe more species and transitions, thus 
to be able to derive reliable abundances and study of the different 
gas components associated with the bow structure. 
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